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We generalise the Buchdahl-Bondi limit for the case of static, spherically symmetric, relativistic 
compact stars immersed in Schwarzschild vacuum in /(A)-theory of gravity, subject to very generic 
regularity, thermodynamic stability and matching conditions. Similar to the case of general relativ¬ 
ity, our result is model independent and remains true for any physically realistic equation of state 
of standard stellar matter. We show that an extra-massive stable star can exist in these theories, 
with surface redshift larger than 2, which is forbidden in general relativity. This result gives a novel 
and interesting observational test for validity or otherwise of general relativity and also provides a 
possible solution to the dark matter problem. 
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I. INTRODUCTION 

In general relativity (GR), differentiable properties and 
regularity conditions of Einstein field equations lead to 
certain interesting bounds on stellar structures [lj. One 
of such fascinating upper bound, is the Buchdahl-Bondi 
bound @0 for the static, spherically symmetric compact 
stars immersed in Schwarzschild vacuum. This bound 
states that the mass to radius ratio 2 M/rt of any regular 
and thermodynamically stable perfect fluid star must be 
less than 8/9. We note that this is a stricter upper bound 
than the Schwarzschild static limit 2 M/rb = 1. In other 
words, even though the mass of the star lies within the 
untrapped region, we cannot have a stable static stellar 
configuration for (8/9)n, < 2 M < rb. The interesting 
points about this general result can be summarised as 
follows: 

(a) This bound can be proved using minimal technical 
conditions, namely the regularity and smoothness of 
metric functions in the stellar interior, matching con¬ 
ditions at the boundary of the star where the space- 
time is smoothly matched to a Schwarzschild exte¬ 
rior. Also for thermodynamic stability, we need to 
impose the condition that the average density is a 
non-increasing function of the radial coordinate. 

(b) This result is model independent. Hence it is true 
for any physically realistic equation of state for the 
stellar matter. 

(c) A direct consequence of this bound is that the grav¬ 
itational redshift z at the stellar surface is bounded 
from above (z < 2) 00. 

Several variations/modifications of this limit have been 
found since, by altering the conditions mentioned above 
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SHllI o r by considering anisotropic stars (see for example 
1214 141] and the references therein). 


Now the key question that arises here is, though this 
result is quite generic within GR, how does it change if 
the theory of gravity is modified. This question is im¬ 
portant even though GR still remains as the successful 
theory of gravity so far, alternative theories of gravity 
do exist. These theories are motivated by the ambiguous 
nature of dark energy in cosmology which is responsible 
for the observed late time accelerated expansion of the 
universe. Another motivation, of course, is the unknown 
nature of dark matter which dominates the matter bud¬ 
get of the galaxies that can be gravitationally detected 
by the galaxy rotation curves. 


To eliminate the need for the so called exotic dark sec¬ 
tors in the universe, an alternative possibility is to con¬ 
jecture that GR is an effective local theory of a more 
general theory on universal scale. Among the numerous 
modifications of GR that naturally provide a late time 
cosmological accelerated expansion, without the need for 
the presence of dark fluids, is f(R)~ gravity. This theory 
is based on a gravitational action that contains an ar¬ 
bitrary but well defined function of the Ricci scalar R 
jl5Hl9f . By expanding the function / around GR, we 
can easily see that the higher order curvature terms nat¬ 
urally admit a phase of accelerated expansion both in 
the early universe as an inflationary phase [20} , and also 
in the late times after passing through a matter domi¬ 
nated decelerating expansion phase[2l|. This theory of 
gravity essentially contain an additional scalar degree of 
freedom which is the Ricci scalar. Hence, this class of 
theories can be considered as a natural extension to GR, 
and contrary to other theories that have the square of 
Ricci or Riemann tensors in the action, the Ostrograd- 
ski instability can be eliminated, despite the equations of 
motion being fourth-order in the metric components. 
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A. The key question 

If we assume the framework of f(R)- gravity as the 
modification to GR on larger scales and strong gravity 
regimes, then the natural question one would ask is: 

Question: Is it possible to have more massive but 
stable compact stars in f (R)-theories of gravity? In 
other words, can the mass to radius ratio of a compact 
star lie in the forbidden region in GR, thus making the 
surface redshift of the star larger than 2? 

The importance of this question is two-fold. If we can 
prove that the mass to radius ratio of a star can be within 
the forbidden region in these theories, then this gives a 
direct experimental test for the validity or otherwise of 
GR in strong gravity regime. Secondly the extra effective 
mass that can be packed in compact stars can account, 
to some extent, for the dark matter in the galaxies, that 
only have gravitational signatures. 


B. In this paper 

In this paper we discuss and try to answer the above 
question, for stable static relativistic compact stars in 
f(R) gravity, immersed in the Schwarzschild vacuum. 
The existence of such stellar solutions is already shown 
in many earlier works (see for example |22l [2?! and the 
references therein). Apart from the regularity, thermo¬ 
dynamic stability and matching conditions that are used 
to prove the Buchdahl-Bondi limit in GR, we now have 
to deal with other extra conditions: 

1. The existence of the Schwarzschild spacetime as a 
solution to the field equations, gives certain condi¬ 
tions on the function /. Since Schwarzschild space- 
time is experimentally well verified around stars 
(e.g. the solar system), we will only consider those 
class of the function / that admits a Schwarzschild 
solution. 

2. For higher order theories such as f(R), we have 
extra matching conditions on the surface of the star 
26, ]29l . [30) , apart from the usual Israel-Darmois 
24 125| conditions in GR. 

3. Since the Ricci scalar is a dynamic degree of free¬ 
dom in this theory, we must impose certain phys¬ 
ically reasonable conditions on it in the interior of 
the star (such as monotonicity) to get a physically 
viable result. 

Taking all these extra conditions into account, in this 
paper we transparently demonstrate the modification of 
Buchdahl-Bondi limit in /(I?)-theories, and show that a 
stable star can exist in the otherwise forbidden region 
in GR, with surface redshift larger than 2. This result 
is absolutely generic since it is model independent as 


in GR. It is true for any equation of state for standard 
stellar matter and also for any function / that satisfies 
the conditions stated above. 

Unless otherwise specified, throughout this paper we 
use units which fix the speed of light and the gravitational 
constant via 8nG = c 4 = 1, and the metric signature is 
+2. 


II. FIELD EQUATIONS FOR RELATIVISTIC 
STATIC STARS IN f(R)- GRAVITY 


To obtain the field equations in /(R)-gravity we begin 
with the modification to the Einstein-Hilbert action by 
generalising the Lagrangian so that the Ricci scalar R 
is replaced by a function f(R). The modified action, 
therefore, is given by 

S= \J dV W~9 f(R) + 2£ M (g a b,ip)\ , ( 1 ) 


where Cm is the Lagrangian density of the matter fields 
ijj, g is the determinant of the metric tensor g a b (a, b = 
0,1, 2, 3), R is the scalar curvature, and f(R) is the real 
valued and well behaved function defining the theory un¬ 
der consideration. Varying the action © with respect to 
the metric over a 4-volume yields the field equations 

Gab f,R - \gab (/ - R f, R ) - VaS7 b f,R + g ab □/,« = T™ , 

( 2 ) 

where f t R = df{R)/dR, □ = V C V C , G a b is the Einstein 
tensor and T is the matter energy momentum tensor 
defined by 


rpM 
1 ab 


2 SC M 

V~9 $g ab 


( 3 ) 


We can easily see that in the special case of f{R) = R , the 
field equations © reduce to the standard Einstein field 
equations. These theories are also known as fourth-order 
gravity, since the term (g a bO—'V a ^b)f,R has fourth-order 
derivatives of the components of the metric tensor. 

We can map the fourth order field equatiions © to the 
effective Einstein equations by 


Gab 




rpM 
-L ab 


rpR _ 

-Lab ~ 


Tab , 


( 4 ) 


where we define T a b as the total energy momentum tensor 
comprising of an effective matter contribution 

fj-iM 

f^ b =-f L , ( 5 ) 

f,R 

and the contribution from the scalar curvature of the 
spacetime as 


rpR _ 

-L ab 


1 

f,R 


-gab (/ — R/,r) + Va^bf,R ~ g a b Of,R 


( 6 ) 
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Writing the field equations in the form Q has an extra 
advantage. We can consider the dynamics of the fourth 
order gravity on the same footing as that of general rel¬ 
ativity with two matter fields, the first is sourced by the 
standard matter and second by the scalar curvature and 
its derivatives. Hence to analyse the properties of these 
field equations, we can implement the well understood 
techniques from general relativity. 


A. Spherical static star in /(P)-gravity 


To study a static spherically symmetric star, immersed 
in the vacuum Schwarzschild exterior, we consider a 
spherically symmetric static metric in the stellar interior 
as 


ds 2 = —c 2 (r)dt , 2 + 


dr 2 


^ 2m[r) j 


r 2 dn 2 , (7) 


where dfl 2 = dd 2 + sin 2 9 d(f> 2 1 c(r) and m{r ) are well de¬ 
fined (at least C 2 ) functions of the radial coordinate. We 
consider the standard matter to be a spherically symmet¬ 
ric perfect fluid with no pressure anisotropy. However we 
cannot impose the same conditions on the “curvature” 
fluid, as in general it is anisotropic. Hence in the above 
coordinate system the total effective energy momentum 
tensor is of the following form: 


Tab = diag \-p{r ), p r (r ), pe (r ), pe (r)] , 

(8) 

where we have set 


. . p M (r) d. . 

P(r) = ± 1 +p R (r) , 

J,R 

(9) 

Pr{r) - +P?(r) , 

J,R 

(10) 

n M ( r ) 

Pe{r) = +Pe(r) ■ 

J,R 

(11) 


Using the definition ([6]) and the metric Q we can easily 
express the pressure anisotropy of the curvature fluid in 
the form 


vf-v 


R 

r 




R'f,RR 



(R 1 f,RR.)' + f ,RR — 
r 


( 12 ) 

where ' denotes differentiation with respect to the coor¬ 
dinate r. 


B. The effective mass 

The important question that arises here is, what 
would be the mass of the star, as measured by a far¬ 


away observer? Since any realistic astrophysical stars 
are immersed in the Schwarzschild vacuum exterior, the 
total effective mass of the star will be equal to the 
Schwarzschild mass in the exterior spacetime. Now from 
the Gg = Xq field equation, we get 

r 

2m' f 

p = —2 “ ^ 2 m(r) = / p(x)x 2 dx . (13) 

o 


Hence we can interpret the function m(r) as the effective 
mass, which is generated by the standard matter and also 
the curvature terms, enclosed within the shell of physical 
radius ‘r’. If r = n, denotes the boundary of the star, 
then m(rb) = M will denote the total effective mass of 
the star, and in the exterior vacuum spacetime this will 
be the Schwarzschild mass. Thus we immediately see 
that the total effective mass of the star in a fourth or¬ 
der theory is different from the mass which is generated 
purely by the standard matter inside the star. The other 
field equations are then used to relate this effective mass 
with the pressure terms. The G\ = X) field equation 
gives 



(14) 


while from the doubly contracted Bianchi identity 
V a G“ = 0 = V a X“, we get 


(cp r y + c'n = — (p e -p r ) ■ (15) 

r 

Substituting equations and (fTT!) in the LHS of (fl5l) , 
we have after simplification 


/ _2m 

IJl - —c' 

= c 

/TON' p B - Pr 

V r dr 

r V r 


Vr 3 / r 

( 


(16) 


C. Problem of an anisotropic star? Not quite. 

The equation (1161) is used extensively in the literature 
to study variations of the Buchdahl-Bondi limit in case of 
stars having pressure anisotropy (see for example El-El 
and the references therein). However this problem (as we 
shall see in this section), is not just another problem of 
an anisotropic star, where the anisotropy parameter is a 
free parameter. In this case, the pressure anisotropy is 
a direct consequence of packaging all the higher deriva¬ 
tive terms in the field equations as the energy momen¬ 
tum tensor of a ‘fictitious’ curvature fluid. Hence, the 
anisotropy parameter depends on the Ricci scalar and 
it’s derivatives, which in turn depend on the metric func¬ 
tion, thereby creating a feedback loop. 

To see this feedback effect more transparently, substi¬ 
tute the curvature terms in (1161) . Since we have already 
considered the standard matter to be a perfect fluid, we 
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can immediately see from equations m and (HD that 
Pe~Pr = pff — p 1 ?- Then using equation (fl2l) in the above 
expression, we get after some simplification 







i 


117) 


The term m/r 3 in the LHS of the above equation denotes 
the average effective density of the star at a radius r. It 
is interesting to note that the metric function c(r) (which 
depends on the gravitational potential) and the function 
f,R play a symmetric role in determining the gradient of 
the effective average density. When f : r> = 1 the above 
relation reduces to the well known relation of general 
relativity jjj, which is used to calculate the Buchdahl- 
Bondi limit. 


III. CONDITIONS ON THE FUNCTIONS 

In order to have a stable and smooth stellar structure 
which can be matched smoothly to the Schwarzschild 
vacuum exterior at the stellar boundary, we have to im¬ 
pose a set of boundary conditions on the interior metric 
functions as well as placing restrictions on the thermo¬ 
dynamic quantities of the star. Furthermore, to avoid 
the presence of ghosts or tachyons and to have the 
Schwarzschild solution as a viable solution of the the¬ 
ory, there are certain conditions on the function f(R). 
In this section we list all these conditions, which will be 
used to calculate the bound on the star mass as dictated 
by the thermodynamic stability. 


A. Conditions of the function f(R) 

To ensure the attractive nature of gravity, that is the 
absence of ghost modes and tachyonic fields we must have 
in the stellar interior 


f,RR> 0. (18) 


B. Matching conditions 


We know that all astrophysical objects are immersed 
in vacuum or almost vacuum spacetime (like any star 
within a stellar system), and hence the exterior spacetime 
around a spherically symmetric star is well described by 
the Schwarzschild geometry. We match two spacetimes 
V ± across the boundary surface E, which in this case will 
be r = rb- The junction surface must be the same in V + 
and V - , which implies continuity of both the metric (the 
first fundamental form) and the extrinsic curvature (th e 
second fundamental form) across E as in GR [U, I25j . 
Moreover, in /(I?)-theories of gravity, continuity of the 
Ricci scalar across the boundary surface and continuity 
of its normal derivative are also required [ID [H, Hi|. 
Matching of the first fundamental form requires that 


c 2 {r b ) 



( 20 ) 


while matching the second fundamental form dictates 
that the total radial pressure at the surface of the star 
must vanish ( p r (r b ) = 0). Therefore using (ThU) we get 


c'{r b ) 


1 M 



( 21 ) 


The extra matching conditions for f(R) gravity give 

R(r b ) = 0 = R'{r b ) . (22) 


C. Regularity Conditions 

The regularity conditions for the interior spacetime 
dictate that all the metric and the thermodynamic func¬ 
tions should be smooth (at least C 2 ) in the interior of the 
star. Hence at the centre of the star, the radial deriva¬ 
tives of all the metric and the thermodynamic functions 
should vanish. Therefore we must have 

c'(0) = p'( 0) = /Lt'(0) = i?'(0) = 0 . (23) 

Also from equation (fl3l) we see that near the centre of 
the star 


m(r) 


(24) 


Furthermore, as proved in [23j , to have the Schwarzschild 
spacetime as a solution of the theory, the function f(R) 
must be at least of class C 3 with 

/(0) — 0 , /,«(0)^0. (19) 

Once these conditions are fulfilled, we can have a star 
with no ghost modes or tachyonic instabilities (which 
are unphysical as they would destroy stable stellar struc¬ 
tures) and furthermore the star can be matched to a vac¬ 
uum Schwarzschild exterior. 


D. Conditions for thermodynamic stability 

From the matching conditions and regularity condi¬ 
tions given above we see that c'(r b ) > 0 and c'(0) = 0. 
Since the function c(r) depends on the gravitational po¬ 
tential, there should be no local extremum of the po¬ 
tential at any non-central shell, for the thermodynamic 
stability of the star. Hence we impose, in the interior 
spacetime, 

c'(r) > 0 , (25) 
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where the equality is only achieved at the central shell. 
Thus the function c(r) will be a monotonically increasing 
function of the coordinate r, ranging from cq = c(0) > 0 
to c{r b ) (which is given by (1201) 1. Furthermore for phys¬ 
ical stability of the stellar structure we also require that 
the total effective density of the star should be a non¬ 
increasing function of r. Therefore we must have 

d /m\ 

Tr A S ° ' < 26 > 

And finally, to avoid the pathologies of f_R = 0 anywhere 
inside the star (which may happen in the Starobinsky 
model and other viable models for negative Ricci scalar), 
we impose a further restriction on the Ricci scalar 

R(r) > 0 . (27) 

Now we can easily see that Rq = -R(O) > 0 while R(r b ) = 
0. For physically viable stellar models we would take the 
function R(r) to be monotonic, which implies 

R'(r) < 0 , (28) 

where the equality is only satisfied at the centre and at 
the surface of the star. 


IV. MODIFICATION OF BUCHDAHL-BONDI 
LIMIT 


In this section, we calculate the bound on the mass 
to radius ratio of a relativistic compact star, subject to 
the stability and regularity conditions (1T51) - (E51) of the 
previous section. As in GR, we shall see that this bound 
is independent of the equation of state of the standard 
matter. Let us denote 


<Kr) = ~J 1- —• (29) 

r V r 

Then we can immediately see that 

= — . (30) 

n V n 

We can now prove easily the following proposition: 

Proposition 1. In the interior spacetime, for any shell 
with the radial coordinate ‘r’, f t n(R(r))c'(r) is positive 
and bounded from below. 

Proof. Using equation (B^l) in ED we get 

f,RT W)c\r)} + c-f [0(r)(/,*)'] < 0 . (31) 

dr dr 

Integrating (by parts) the above equation from r to r b , 
where r b is the boundary of the star, we obtain 

[f,RW)c% b + [cm(f, R y]] r ; 


rb 

-2 J f, rrR' c'cj){r) < 0 . 


r 


Since in the interior spacetime R' < 0, the third term 
in LHS of the above equation has a positive contribution 
and hence dropping that term does not change (rather 
strengthens) the inequality. Therefore we have 

f,R{R(n)) [<Kn)c'(r 6 )] - f, R (R(r)) [0(r)c'(r-)] 

+c(r b ) [(t>{r b )f^ rrR'{ r b )\ - c(r) [<j)(r) f, rrR'] < 0. 

(33) 

Now using (1221) in the above equation, the third term 
vanishes, while by C51) the fourth term has a positive 
contribution to the LHS and can be dropped without 
altering the inequality. Substituting c'(r b ) from (l2lT) we 
get 


f,R(R(r))£ > 
dr 


fA 0) Mr 
rl 


(34) 

□ 


Since the above inequality is true Vr e [0,rt,], we can 
integrate the above from the centre to the surface, with¬ 
out changing the inequality. Integrating we get 

rb 

f,R(0)c(r b ) - f, r(Ro)cq - J f.R R R'cdr > 

o 


fA0)M f 

rl J 

o 



(35) 


To get a bound on the integral term above, we state and 
prove the following proposition: 

Proposition 2. For any thermodynamically stable stel¬ 
lar model, satisfying the conditions \1S\ )~ \28\) . the integral 
of (l/cj)(r)) from the centre to the surface is bounded from 
below, and the bound is equal to (rj)/2M) [1 — r b (f{r b )] . 

Proof. To prove this, we follow the same steps as given 
in [7|. Let us write the total effective density function 
/r(r) as a sum of constant average density no = 6M/rl 
and p{r) as the variation on the average, that is 

. . 6 M 

hi r ) = — + P( r ) ■ (36) 

r b 


This implies 


r 

r 3 I 

2 m(r) = 2 M— + / p(r)r 2 dr . (37) 

r h J 


From the above equation we can easily see that 


rb 

J p{r)r 2 dr = 0 . 


(32) 


o 


(38) 
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Also since the effective density is a monotonically non 
increasing function from the centre to the surface of the 
star, we must have the central density greater than the 
average density. Therefore we have p(0) > 0 and p' < 0. 
The integral term in (l37l) is always positive and goes to 
zero as r —> r b . This implies, Vr £ [0, r&] that 


and therefore 


2 M^t < 2 to , 

n 


(39) 


/ 




s A. 

2 M 


[1 - r b <j>(r b )\. 


(40) 


□ 


We can use the above equation in the RHS of (1351) . 
without altering the inequality. Hence we now get 


rb 

/,fl(0)c(rb) - f,n(Ro)co - J f, RR R'cdr > 

o 


fA Q) 

2 



(41) 


Proof. From the regularity conditions at the centre of the 
star, and also by the conditions on the function f(R), 
we can see that the RHS of the equation (PHI) is strictly 
greater than zero. Therefore we get 


f, 2M r/,*(0) 

V 1_ Th“ 


fARo) 


fA°) 

2 


>o, 


(45) 


which, by squaring the LHS, can be simplified to the 
following expression 


i , n fA R o )] 2 

+ /.*(0) J 



> 1 . 


(46) 


Taking the quantity 2 M to the RHS and simplifying we 
finally obtain the upper bound as 


2 M < 


A f,R{R o) 

4 fM 0) 


1 + 


fMRo ) 


fMo) 


1 + 2 


fMRo) 


fM o) 


r b 


(47) 


By our assumptions on the Ricci Scalar and the function 
f(R), we can immediately see that > 1- When 

this factor equals unity, we regain the usual Buchdahl- 
Bondi upper bound 2 M < (8/9 )r b . When this term is 
much larger than unity, this upper bound tends to the 
Schwarzschild static limit 2 M = r b - □ 


Now consider the integral term in the LHS of the above 
equation. As R' < 0, we have 


- / f,RRR'c{r)d = 


f,RRR'c{r)dr 


(42) 


Since this term has a positive contribution in the LHS 
of (1411) . we cannot drop it. However as we have seen 
c'(r) > 0 throughout the interior of the star, therefore 
c(r b ) > c(r), Vr £ [0, r b ] and hence the inequality will not 
change if we substitute 


The above theorem ensures that if we suitably modify 
general relativity by a function f{R) ^ R, we can have 
a thermodynamically stable, regular spherical star with 
an effective mass M in the region (8/9)r;, < 2 M < r b , 
which is the forbidden region in general relativity. Hence 
in these theories we can have more massive but stable 
compact stars, that may be one of the keys to the dark 
matter problem. Also for these stars the surface redshift 
z > 2, which can in principle be detected experimentally. 
To see this effect more transparently, let us assume a 
small modification from general relativity. In other words 
let 


c(r b ) 


J.rrR' dr 


= c(r b ) (fAR o) - /+(0)), (43) 


o 

for the integral term in (I4T1) . On substitution, and using 
we finally have 


1 - 


2 M 
r b 


fA°) 


I,r{Ro) 


fA°) 

2 


> fARoA ■ 


(44) 

We are now in a position to state and prove the following 
theorem: 


Theorem 1. The regularity and thermodynamic stabil¬ 
ity conditions on a spherically symmetric and static star 
of radius r b , immersed in the Schwarzschild vacuum in 
f(R)-gravity, impose an upper bound on the total effec¬ 
tive mass of the star. This upper bound lies between 
the Buchdahl-Bondi bound for general relativity and the 
Schwarzschild static limit 2 M = r b . 


fARd) 

fA°) 


— 1 + Of. , 


0 < a « 1 . 


(48) 


Then the upper bound on the effective mass of the star 
becomes 


2M < U i+ ^) n - (49) 

Hence the net fractional increase of the effective mass, 
defined as SM = ( M e ff — Mgr)/Mgr (where Mgr is 
the upper bound in general relativity), can be written as 


SM _ 4a 
Mgr 27 


(50) 


This may easily and naturally account for the extra mas¬ 
sive neutron stars in the sky. We can immediately cal¬ 
culate the maximum surface redshift for such stars. We 
get 


2 + Zmax — 2(1 + a) 


(51) 
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Hence, from our observational data from the compact 
stars, we can in principle experimentally verify any devi¬ 
ation from general relativity in the high curvature regime. 

V. DISCUSSION 

In this paper, we studied model independent bounds on 
spherically symmetric stellar structures in /(I?)-gravity. 
In particular, our results are independent of the matter 
distribution and the equation of state. Subject to very 
generic conditions of regularity and thermodynamic sta¬ 
bility in the interior of the star and the matching condi¬ 
tions at the stellar surface to the vacuum Schwarzschild 
exterior, we transparently demonstrated that the mass to 
radius ratio of the star is bounded from above, and this 
is a stricter bound than the Schwarzschild static limit. In 
other words, we generalised the Buchdahl-Bondi bound 
on static stars in GR to f(R)- theories. 

We also showed that this upper bound is larger than 
the Buchdahl-Bondi limit of GR, whenever f(R ) ^ R. 
Hence, in principle, we can pack extra effective mass in a 
stable compact star in these theories, which is forbidden 
in GR. These extra massive stars may be one of the solu¬ 
tions to the dark matter problem, that manifests through 


the rotation curves of the galaxies. 

Furthermore, as a direct consequence of this extra ef¬ 
fective mass, we proved that the surface redshift of an 
extra-massive compact star can be greater than 2, which 
is the upper limit in the case of general relativity. This 
gives a novel and interesting scenario, where we can ob- 
servationally verify the validity or otherwise of general 
relativity in the strong gravity/ high curvature regime. 

An interesting possibility of further study in this sce¬ 
nario would be the violation of Chandrasekhar limit. We 
know, it is possible that Chandrasekhar limit may be vio¬ 
lated in the presence of strong magnetic fields. The exis¬ 
tence of super-Chandrasekhar white dwarfs has been mo¬ 
tivated on the grounds of a polytropic equation of state in 
GR for an anisotropic matter distribution a It would 
be interesting to see, how the curvature terms in f{R)~ 
gravity modify this limit. 
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